Heat transfer between elastic solids with randomly rough surfaces 
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We study the heat transfer between elastic solids with randomly rough surfaces. We include both 
the heat transfer from the area of real contact, and the heat transfer between the surfaces in the non- 
contact regions. We apply a recently developed contact mechanics theory, which accounts for the 
hierarchical nature of the contact between sohds with roughness on many different length scales. For 
elastic contact, at the highest (atomic) resolution the area of real contact typically consists of atomic 
(nanometer) sized regions, and we discuss the implications of this for the heat transfer. For solids 
with very smooth surfaces, as is typical in many modern engineering applications, the interfacial 
separation in the non-contact regions will be very small, and for this case we show the importance of 
the radiative heat transfer associated with the evanescent electromagnetic waves which exist outside 
of all bodies. 



1. Introduction 

The heat transfer between solids is a topic of great im- 
portance. Classical applications include topics such as 
cooling of microelectronic devices, spacecraft structures, 
satellite bolted joints, nuclear engineering, ball bearings, 
tires and heat exchangers. Other potential applications 
involve microelectromechanical systems (MEMS). Heat 
transfer is also of crucial importance in friction and wear 
processes, e.g., rubber friction on hard and rough sub- 
strates depends crucially on the temperature increase in 
the rubber-countersurface asperity contact regions 

A large number of papers have been published on 
the heat transfer between randomly rough surfaces 
However, most of these studies are based on asperity 
contact models such as the model of Greenwood and 
WiUiamson (GW)[i]. Recent studies have shown that 
the GW- model (and other asperity contact models [1|) 
are very inaccurate 0, mainly because of the neglect 
of the long-range elastic coupling^. That is, if an asper- 
ity is pushed downwards somewhere, the elastic deforma- 
tion field extends a long distance away from the asperity, 
which will influence the contact involving other asperi- 
ties further awayQ. This effect is neglected in the GW 
theory, but it is included in the contact mechanics model 
of Perssonfol. [loL H, 12, l^], which we use in the present 
study. In addition, in the GW model the asperity contact 
regions are assumed to be circular (or elliptical) while the 
actual contact regions (at high enough resolution) have 
fractal-like boundary hnes[ll, Q, [S], see Fig. [T] Thus, 
because of their complex nature, one should try to avoid 
to directly involve the nature of the contact regions when 
studying contact mechanics problems, such as the heat or 
electric contact resistance. The approach we use in this 
paper does not directly involve the nature of the contact 
regions. Finally, we note that for elastically hard solids 
the area of real (atomic) contact A may be a very small 
fraction of the nominal or apparent contact area Aq, even 
at high nominal squeezing pressures H, 17|. 

Another important discovery in recent contact me- 
chanics studies is that for elastic contact, the contact 




FIG. 1: The black area is the contact between two elastic 
solids with randomly rough surfaces as obtained using molec- 
ular dynamics. For surfaces which have fractal-like roughness 
the whole way down to the atomic length scale, the contact 
at the highest magnification (atomic resolution) typically con- 
sists of nanometer-sized atomic clusters. Adapted from Ref. 

El. 



regions observed at atomic resolution may be just a few 
atoms wide, i.e., the diameter of the contact regions may 
be of the order of ~ 1 nmfU 3 2^- The heat transfer 
via such small junctions may be very different from the 
heat transfer through macroscopic sized contact regions, 
where the heat transfer usually is assumed to be propor- 
tional to the linear size of the contact regions (this is also 
the prediction of the macroscopic heat diffusion equa- 
tion), rather than the contact area. In particular, if the 
typical phonon wavelength involved in the heat transfer 
becomes larger than the linear size of the contact regions 
(which will always happen at low enough temperature) 



the effective heat transfer may be strongly reduced. Sim- 
ilarly, if the phonons mean free path is longer than the 
linear size of the contact regions, ballistic (phonon) en- 
ergy transfer may occur which cannot be described by 
the macroscopic heat diffusion equation. These effects 
are likely to be of crucial importance in many modern 
applications involving micro (or nano) sized objects, such 
as MEMS, where just a few atomic-sized contact regions 
may occur. However, for macroscopic solids the ther- 
mal (and electrical) contact resistance is usually very in- 
sensitive to the nature of the contact regions observed 
at the highest magnification, corresponding to atomistic 
(or nanoscale) length scales. In fact, the heat transfer 
is determined mainly by the nature of the contact re- 
gions observed at lower magnification where the contact 
regions appear larger (see Sec. 5 and 21, 2§|), see Fig. 
[2l For example, in Sec. 2.2.1 we show that for self- 
affine fractal surfaces the contact resistance depends on 
the range of surface roughness included in the analysis as 
~ r{H) — {qo/qi)^ , where go and qi are the smallest and 
the largest wavevector of the surface roughness included 
in the analysis, respectively, and H is the Hurst exponent 
related to the fractal dimension via D[ = 3 — H. The 
number r{H) depends on H but is of the order of unity. 
In a typical case H « 0.8, and including surface rough- 
ness over one wavevector decade qo < q < qi = lOqo 
results in a heat resistance which typically is only ~ 10% 
smaller than obtained when including infinitely many 
decades of length scales (i.e., with qi = oo x qo). At 
the same time the area of real contact approaches zero 
as qo/qi 0. Thus, there is in general no relation be- 
tween the area of real contact (which is observed at the 
highest magnification, and which determines, e.g., the 
friction force in most cases), and the heat (or electrical) 
contact resistance between the solids. One aspect of this 
in the context of electric conduction was pointed out a 
long time agofisj: if an insulating film covers the solids 
in the area of real contact, and if electrical contact occurs 
by a large number of small breaks in the film, the resis- 
tance may be almost as low as with no film. Similarly, 
the thermal contact resistance of macroscopic solids usu- 
ally does not depend on whether the heat transfer occur 
by diffusive or ballistic phonon propagation, but rather 
the contact resistance is usually determined mainly by 
the nature of the contact regions observed at relative low 
magnification. 

Note that as H decreases towards zero (or the frac- 
tal dimension D[ 3) one needs to include more and 
more decades in the length scales in order to obtain the 
correct (or converged) contact resistance, and for H = 
(or D[ = 3) it is necessary to include the roughness on 
the whole way down to the atomic length scale (assum- 
ing that the surfaces remain fractal-like with H — the 
whole way down to the atomic length scale). Most nat- 
ural surfaces and surfaces of engineering interest have (if 
self-affine fractal) H > 0.5 (or Df < 2.5), e.g., surfaces 





FIG. 2: The contact region (black area) between two elas- 
tic solids observed at low (left) and high (right) magnifica- 
tion. The contact resistance depends mainly on the long- 
wavelength roughness, and can usually be calculated accu- 
rately from the nature of the contact observed at low magni- 
fication (left). 



prepared by crack propagation or sand blasting typically 
have H « 0.8, and in these cases the contact resistance 
can be calculated accurately from the (apparent) con- 
tact observed at relatively low magnification. However, 
some surfaces may have smaller Hurst exponents. One 
interesting case is surfaces (of glassy solids) with frozen 



capillary waves [ij, \2M (which are of great engineering 
importance [241 )■ which have H = 0. The heat transfer 



between such surfaces may be understood only by study- 
ing the system at the highest magnification correspond- 
ing to atomic resolution. 

In this paper we will consider the heat transfer between 
(macroscopic-sized) solids in the light of recent advances 
in contact mechanics. We will study the contribution to 
the heat transfer not just from the area of real contact 
(observed at atomic resolution), but also the heat transfer 
across the area of non-contact, in particular the contri- 
bution from the fluctuatin g el ectromagnetic field, which 
surrounds all solid obiectsjisl. [2^ . For high-resistivity 
materials and for hard and very flat surfaces, such as 
those involved in many modern applications, e.g., MEMS 
applications, this non-contact radiative heat transfer may 
in fact dominate in the total heat transfer (at least un- 
der vacuum condition). We note that for flat surfaces (in 
vacuum) separated by a distance d larger than the ther- 
mal length (It — cfi/k^T^ the non-contact heat transfer 
is given by the classical Stefan-Boltzman law, and is inde- 
pendent of d. However, for very short distances the con- 
tribution from the evanescent electromagnetic waves to 
the heat transfer will be many orders of magnitude larger 
than the contribution from propagating electromagnetic 
waves (as given by the Stefan-Boltzman law) [2^. 

In most applications (but not in spacecraft applica- 
tions) one is interested in the heat transfer between solid 
objects located in the normal atmosphere and sometimes 
in a fluid. Most solid objects in the normal atmosphere 
have organic and water contamination layers, which may 
influence the heat transfer for at least two reasons: (a) 
Thin (nanometer) contamination layers may occur at the 



interface in the asperity contact regions, which will ef- 
fect the acoustic impedance of the contact junctions, and 
hence the propagation of phonon's between the solids 
(which usually is the origin of the heat transfer, at least 
for most non-metallic systems), (b) In addition, cap- 
illary bridges may form in the asperity contact regions 
and effectively increase the size of the contact regions 
and increase the heat transfer. In the normal atmosphere 
heat can also be transferred between the non-contact re- 
gions via heat diffusion or (at short separation) ballistic 
processes in the surrounding gas. For larger separations 
convective processes may also be important. 

In the discussion above we have assumed that the solids 
deform elastically and we have neglected the adhesional 
interaction between the solids. The contact mechanics 
theory of Persson can also be applied to cases where ad- 
hesion and plastic flow are important, and we will briefly 
study how this may affect the heat transfer. Most solids 
have modified surface properties, e.g., metals are usually 
covered by thin oxide layers with very different conduc- 
tivities than the underlying bulk materials. However, as 
mentioned above, this may not have any major influence 
on the contact resistance. 

Recently, intense research has focused on heat transfer 
through atomic or molecular-sized junctions [27l. [28j . In 
light of the discussion presented above, this topic may 
also be important for the heat transfer between solids, 
because of the nanometer-sized nature of the contact re- 
gions between solids with random roughness. 

This paper is organized as follows: In Sec. 2 we de- 
scribe the theory for heat transfer between two solids 
with randomly rough surfaces. We consider both the 
heat flow in the area of real contact, and between the 
surfaces in the non-contact area. Sec. 3 presents a short 
review of the contact mechanics theory which is used to 
obtain the quantities (related to the surface roughness) 
which determine the heat transfer coefficient. In Sec. 4 
we present numerical results. In Sec. 5 we discuss the 
influence of plastic flow and adhesion on the heat trans- 
fer. Sec. 6 presents an application to the heat transfer 
between tires and the air and road surface. In Sec. 7 we 
discuss a new experiment. In Sec. 8 we present experi- 
mental results. In Sec. 9 we point out that the developed 
theory can also be applied to the electric contact resis- 
tance. Sec. 10 contains the summary and conclusion. 
Appendix A-E present details related to the theory de- 
velopment and some other general information relevant 
to the present study. 

2. Theory 

2.1 Heat transfer coefficient 

Consider two elastic solids (rectangular blocks) with 
randomly rough surfaces squeezed in contact as illus- 
trated in Fig. [3l Assume that the temperature at the 
outer surfaces z = —do and z — di is kept fixed at Tq 
and Ti , respectively, with Tq > Ti. Close to the interface 
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FIG. 3: Two elastic solids with nominally flat surfaces 
squeezed together with the nominal pressure po. The heat 
current Jz(x) at the contacting interface varies strongly with 
the coordinate x = {x, y) in the xy-p\&ne. The average heat 
current is denoted by Jo = (Jz(x)). 



the heat current will vary rapidly in space, J = J(x, z), 
where x = {x,y) denote the lateral coordinate in the xy- 
plane. Far from the interface we will assume that the heat 
current is constant and in the z-direction, i.e., J = Jqz. 
We denote the average distance between the macro as- 
perity contact regions by A (see Ref. [l^). We assume 
that A << L, where L is the linear size of the apparent 
contact between the elastic blocks. The temperature a 
distance ~ A from the contacting interface will be approx- 
imately independent of the lateral coordinate x — (x, y) 
and we denote this temperature by Tg and T[ for z = — A 
and z — \, respectively. The heat current for |z| >> A is 
independent of x and can be written as (to zero order in 
A/do and A/di): 



^0 



do 



-Kl- 



T i-Tj 
di 



(1) 



where kq and ki are the heat conductivities of the two 
solid blocks. We assume that the heat transfer across the 
interface is proportional to Tq — T[ and we define the heat 
transfer coefficient a so that 



Jo = a(T^ - Ti) 
Combining (1) and (2) gives 

To 



Jo = 



do Ho 



(2) 



(3) 



This equation is valid as long as A << L and A << 
do, di. Note that a depends on the macroscopic (or 
nominal) pressure which act at the interface. Thus if the 



macroscopic pressure is non-uniform, as is the case in 
many practical applications, e.g., when a ball is squeezed 
against a flat, one need to include the dependence of a 
on X. Thus in general 



J(x)=a(x) [T^(x)-T;(x)] 



(4) 



One expect the contribution to a from the area of real 
contact to be proportional to the heat conductivity k 
(for simplicity we assume here two solids of the same 
material). Assuming only elastic deformation, contact 
mechanics theories show that for low enough squeezing 
pressure po, the area of real contact is proportional to po; 
and the size distribution of contact regions (and the in- 
terfacial stress probability distribution) are independent 
of pq. Thus one expect that a is proportional to po- For 
randomly rough surfaces the contact mechanics depends 
only on the (effective) elastic modulus E* and on the sur- 
face roughness power spectrum C (g) . Thus the only way 
to construct a quantity which is proportional to pqk, and 
with the same dimension as Jq/ AT, using the quantities 
which characterize the problem, is 

PqK 



E*uo 



where Uq is a length parameter which is determined from 
the surface roughness power spectrum C{q). For self- 
affine fractal surfaces, C{q) depends only on the root- 
mean-square roughness /inns, the fractal dimension D{ 
which is dimension less, and on the low and high cut- 
off wavevectors qo and qi. Thus in this case uq — 
Qohrias)- This result is consistent with 
the analysis presented in Sec. 2.2.1. Using the GW- 
theory result in an expression for a of the form given 
above, but with a different function / which now (even 
for low squeezing pressures) also depends on po/E* (see, 
e.g., Ref. |29|]). 

2.2 Calculation of a 

The heat current J and the heat energy density Q are 
assumed to be given by 



J = -kVT, 



Q = pCvT 



where k is the heat conductivity, p the mass density and 
Cv the heat capacitivity. We consider a steady state 
condition where Q is time independent. Thus the heat 
energy continuity equation 



reduces to 







We assume that the surface roughness at the interface 
is so small that when solving the heat flow equation we 
can consider the surfaces as flat. However the heat flow 



across the interface will be highly non-uniform and given 
by the heat current Jz(x) (we assume |V/i| << 1, where 
/i(x) is the surface height profile). Let us first study the 
heat flow in the upper solid. We can take into account 
the heat flow from the lower solid by introducing a heat 
source at the interface z = i.e. 



V^T = -2Jz(x)(5(z)/ki 



(5) 



Similarly, when studying the temperature in the lower 
solid we introduce a heat sink on the surface z — so 
that 

V2t = 2Jz(x)(5(z)/ko (6) 
Let us first study the temperature for z > 0. We write 



J4x) = J d\ JMV^'^ 

JM) = J cP^ J.(x)e-'i- 

From (5) we get 



(7) 



(8) 



T(x, z) =Ti Jo(z - di) 



x-\-kz) 



(9) 



TTKl J —(p- — 

where Jq = {Jz{Ti)) is the average heat current and 

AJ,(x) = J,(x)- Jo (10) 
Performing the fc-integral in (9) gives 

T(x, z) = Ti - —Jo{z - di) 



-— [ d^q -AJ,{q)e''i-''-'"' 
Kl J q 



(11) 



Similarly, one obtain for the temperature field for z < 0: 
T(x, z) = To - —Jo{z + do) 

Kq 



-— [ d^q -A.h{(i)e''i ''+'" 

Ko J q 



fl2) 



Let us define 

V'(x) =T(x,-0)-T(x,+0) 
Using (11) and (12) we get 



[d'q -AUq)e'- 

kJ q 



where 



11 1 

K Kq Kl 



From (13) we get 



V(q)-A«(q)-— AJ,(q) 
Kq 



where 



M = To-T,-(^ + ^)jo 
We will now consider two different cases: 



(13) 



(14) 



(15) 



(16) 



2.2.1 Heat flow through the area of real contact 

Let us consider the area of real contact. In the contact 
region Jz(x) will be non-zero but = T(x, +0) — 

r(x, — 0) will vanish. On the other surface area Jz{^) 
will vanish. Thus we must have 

J,(x)V(x) =0 

everywhere. This implies 

j d\' JM~c{)^{c() = (17) 

for all q. Combining (15) and (17) gives 

MJ,(q) ~ - [ d\'^JM-^)^JM) = 
K J q' 

The ensemble average of this equation gives 

M{JM)) -\j d'g'^(J.(q- q')AJ.(q')) = (18) 

From (8) we get 

(J,(q = 0)) = (27r)-2AoJo. 
Thus the q = component of (18) gives 

MAoJo - ^ / d\-{\^JM)?) = (19) 
K J q 

where Aq is the nominal contact area. Combining (16) 
and (19) and solving for Jq gives an equation of the form 
(3) with 



1 - (27^)' 1 
a K AqJq 



d\-^{\^JM)?) (20) 



We now assume that the heat current at the interface 
is proportional to the normal stress: 



We can also write (21) as 



Jz(x)/Jo « CTz(x)/po, 



(22) 



where po is the average pressure. We note that (22) im- 
plies that the current density Jz(x) will be non- vanishing 
exactly where the normal stress ^^(x) is non- vanishing, 
which must be obeyed in the present case, where all the 
heat current flow through the area of real contact. We 
note that the heat transfer coefficient depends mainly 
on the spatial distribution of the contact area and this 
is exactly the same for the pressure distribution cr(x) as 
for the current distribution Jz(x). Thus the fact that in 
a particular asperity contact region the pressure cr(x) is 
not proportional to Jz{^) is not very important in the 
present context (see Appendix A and below). 
Substituting (22) in (20) gives 



1 



a 



K AqP^ 



\ [ d^q-{\Aaz{ci) 



We can write 



where 



d\-{\Aaiq)\' 



(23) 



(24) 



(25) 



is the stored elastic energy per unit (nominal) surface 
area[l3|- In (25) E* is the effective elastic modulus 



1 

E* 



En 



E, 



where Eq and Vq are the Young's elastic modulus and 
the Poisson ratio, respectively, for solid and similar for 
solid 1. We have shown elsewhere that for small enough 
load[l3] Ue\ ~ uopa where uq is a length of order the 
root-mean-square surface roughness amplitude. Thus 



Pqk 
E*uo' 



(26a) 



Note that for small load the squeezing pressure po de- 
pends on the (average) interfacial separation u via the 
exponential law po exp(— u/uo). Thus the vertical 
stiffness dpo /du — —po / uq so we can also write 



K dpo 
'^~dU' 



(266) 



Jz(x) « H<7z{x.) 



(21) 



This equation is, in fact, exact (see Appendix B and Ref. 
p^). which shows that the heat transfer is mainly deter- 
mined by the geometrical distribution of the contact area 
(given by the region where azi'x) is non- vanishing) , and 
by the thermal interaction between the heat flow through 
the various contact spots (see Appendix A). 



The length parameter uq in (26a) can be calculated 
(approximately) from the surface roughness power spec- 
trum C{q) using [m 



where 



Mo 



w(q) 



Vtt / dq q^C{q)w{q) 



q X -1/2 

TT / dq'q'^C{q') 
90 



where go is the long-distance cut-off (or roll-off) wavevec- 
tor and qi the wavevector of the shortest wavelength 
roughness included in the analysis. Assume that the com- 
bined surface roughness is self afhne fractal for qo < q < 
qi. In this case 



C{q) = 



2(H+1) 



where H is the Hurst exponent related to the fractal di- 
mension via D{ — 3 ~ H . Substituting this C{q) into the 
equations above gives 



Uo ! 

where 
r{H) 



2(1 -g) 
nH 



H 

2(1 -i^) Ji 



1/2 



r{H) - ( 

91 



H 



dx (a;-l)-^/'x-i/P(i-^)l 



Note that r{H) is of order unity (see Ref. \l(J\)- As dis- 
cussed in the introduction this implies that the contact 
resistance in general is determined accurately by one or 
two decades of the longest-wavelength roughness com- 
ponents, and that there is no relation between the area 
of real contact (which is observed at the highest mag- 
nification, and which determines, e.g., the friction force 
in most cases), and the contact resistance between the 
solids. 

Note that from (3) it follows that one can neglect the 
heat contact resistance if 

K/d << a 

where n/d is the smallest of Ko/do and Ki/di. Using (25) 
this gives 

d » uo{E*/pq) 

We note that in modern high-tech applications the linear 
size (or thickness) d of the physical system may be very 
small, and in these cases the contact heat resistance may 
be particular important. 

If roughness occurs only on one length scale, say with 
wavelength A and height h, then the pressure necessary 
for complete contact will be of order 

po « E*hl\ 



Substituting this in (26a) gives 



(27) 



where we have used that uq sa h. Thus, a"^ ~ \k~^ 
which is the expected result because the denominator in 
(3) is only accurate to zero order in \k~^ . [Alternatively, 
substituting (27) in (3) gives a term of the type {d+X)n^'^ 
which is the correct result since d in (3) should really be 
A.] 

As an example [30t. consider two nominal flat steel 
plates (in vacuum) with the thickness — di = 0.5 cm 
and with the root-mean-square roughness ~ 1 /xm. The 
plates are squeezed together with the nominal pressure 
Po = 0.1 MPa. The ratio between the measured surface 
and bulk thermal contact resistance is about 150. Using 
(3) we get 

AT/ Jo = 2doKo^ 

Thus, the (theoretical) ratio between the surface and the 
bulk contributions to the thermal resistance is: 

Ko 



2adQ ' 



where kq is the heat conductivity of the bulk steel. Using 
(25) with K = Ko/2 this gives 



2adn 



uoE*_ 
do Pa 



(28) 



With (from theory) wq w 1 /^m, and E* w 110 GPa, po = 
0.1 MPa and 2do = 1 cm, from (28) the ratio between 
the thermal surface and bulk resistance is w 200, in good 
agreement with the experimental data. 

The discussion above assumes purely elastic deforma- 
tions. However, plastic flow is likely to occur in the 
present application at short enough length-scales, ob- 
served at high magnification. Since the heat flow is deter- 
mined mainly by the long-wavelength roughness compo- 
nents, i.e., by the roughness observed at relative low mag- 
nification, when calculating the heat transfer one may 
often assume that the surfaces deform purely elastically, 
even if plastic deformation is observed at high magnifi- 
cation, see Sec. 5. 

2.2.2 Heat flow through the non-contact area 

Let us now assume that 

J,(x) = /3(x) [T(x, -0) - r(x, -fO)] = /3(x)V'(x) 
From (15) we get 

^(q) = M5(q) 



/ d^q' /3(q-q') 

nq J 



^0 



^(q') (29) 



Next, note that 



Jo = ^ d^x J,(x) = ^^j<^x /3(x)V'(x) 



d\ /3(-q)^(q) 



(30) 



Eq. (29) can be solved by iteration. The zero-order 
solution 

V;(q) = M5(q) 
Substituting this in (30) gives 

Jo = M^^/3(q = 0) = Mp (31) 

where 

= {Pi^)) = l-Jd'x /3(x) 

is the average of /3(x) over the whole interfacial area Aq. 
Substituting (16) in (31) and solving for Jq gives an equa- 
tion of the form (3) with a = f3. 
The first-order solution to (29) is 



V(q) = M5(q) - — /3(q) 



1 



^0 



^(q) 



(32) 



Substituting (32) in (30) gives again an equation of the 
form (3) with 



a = j3- 



kAo 



J d2gi(|/3(q)H 



1 



5(q) 



, (33) 



where we have added (..) which denotes ensemble aver- 
age, and where we used that 

(/3(q)/3(-q))=(|/3(q)|2) 

We can rewrite (33) as follows. Let us define the corre- 
lation function 



1 



d^x (/3(x)/3(0))e''»-= 



Note that 



C,(q) = ^(|/3(q)l^) 



(34) 



(35) 



This equation follows from the fact that the statistical 

properties arc assumed to be translational invariant in 
the x-plane, and is proved as follows: 

^^^"^^ = (2^)^ /^'"^ (/3(x)/3(0))e^^- 



(2 



^ j d^x (/3(x + x')/3(x'))e^'>- 



= / d-'x" (/3(x")/3(x'))e^^-(-"--') 

This equation must be independent of x' and we can 
therefore integrate over the x'-plane and divide by the 
area Aq giving 

CM) = / d^x'd^x" {(3{^W))e'^<-"-') 



Aq 



(l/3(q)l') 



Let us define 



We get 



A/3(x) = ;3(x) - p 
A/3(q) = /3(q) - /3,5(q) 



(36) 



and thus 

(|A/?(q)|^) = 
where we have used that 



^0 



(37) 



^<5(q) = i^/3(q)5(q) 

and that 

W(-q) = '5(q)(^ / d'x e-^^ - 
Using (33) and (37) gives 

a = p-- I d\q-'CAM) 

K J 

Let us write 

(A;3(x)A/3(0)) = ((A/3)2)/(x) 
where /(O) = 1. We write 

/(x) = J d\ /(q)e^^-- 

so that /(x = 0) = 1 gives 

j d\ fid) = 1 

Using (39) and (40), Eq. (38) takes the form 

a = P - {{Apf)K-H 

where the correlation length 

,_ /Ag-V(q) 
Jd^qfici) 
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FIG. 4: Solid line: The calculated [using (42)] heat current 
per unit area, Jo, between two (amorphous) sihcon dioxide 
bodies, as a function of the temperature difference AT. The 
solids have flat surfaces separated by d = 1 nm. One solid 
is at the temperature T = 296 K and the other at T + AT. 
Dashed line: linear function with the slope given by the initial 
slope (at AT = 0) of the solid line. 



For randomly rough surfaces with isotropic statistical 
properties /(q) depends only on q = |q| so that 



I 



/o°° dq /(g) 
/o°° dq qf{q) 



Most surfaces of engineering interest are fractal- like, with 
the surface roughness power spectrum having a (long- 
distance) roll-off wavevector go- In this case one can show 
from that I q^^. For the surface used in the numerical 
study presented below in Sec. 4 one have qo ~ 10^ m^^ 
(see Fig. [8]). Furthermore, in this case (for amorphous 
silicon dioxide solids) k « 1 W/mK and if we assume 
that ((A/3)^) is of order we get the ratio between the 
second and the first term in (41) to be of order /?/ (qon) « 
0.01, where we have used that typically (see Fig. [5]) /3 « 
0.1 MW/m^K. Thus, in the application presented in Sec. 
4 the second term in the expansion (41) is negligible. 

Eq. (41) represent the first two terms in an infinite 
series which would result if (29) is iterated to infinite 
order. The result (41) is only useful if the first term /3 
is much larger that the second term. If this is not the 
case one would need to include also higher order terms 
(in principle, to infinite order) which becomes very hard 
to calculate using the iterative procedure. By comparing 
the magnitude between the two terms in (41) one can 
determine if it is legitimate to include only the lowest 
order term (3. 

We now consider two applications of (41), namely the 
contribution to the heat transfer from (a) the electro- 
magnetic field (in vacuum) and (b) from heat transfer 
via a gas (e.g., the normal atmosphere) which we assume 
is surrounding the two solids. 

(a) Radiative contribution to a (in vacuum) 



The heat flux per unit area between two black-bodies 
separated hy d » dx = ch/kBT is given by the Stefan- 
Boltzmann law 



60h 



where Tq and Ti are the temperatures of solids 1 and 
2, respectively, and c the light velocity. In this limiting 
case the heat transfer between the bodies is determined 
by the propagating electromagnetic waves radiated by 
the bodies and does not depend on the separation d be- 
tween the bodies. Electromagnetic waves (or photons) 
always exist outside any body due to thermal or quan- 
tum fluctuations of the current density inside the body. 
The electromagnetic field created by the fluctuating cur- 
rent density exists also in the form of evanescent waves, 
which are damped exponentially with the distance away 
from the surface of the body. For an isolated body, the 
evanescent waves do not give a contribution to the energy 
radiation. However, for two solids separated hy d < dx, 
the heat transfer may increase by many orders of mag- 
nitude due to the evanescent electromagnetic waves-this 
is often referred to as photon tunneling. 

For short separation between two solids with flat sur- 
faces {d << dx), the heat current due to the evanescent 
electromagnetic waves is given bvpsj 



Jn = 
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(42) 



where 



and 



e(t^) + 1 



where e(w) is the dielectric function. From (42) it follows 
that the heat current scale as l/d"^ with the separation 
between the solid surfaces. The heat current is especially 
large in the case of resonant photon tunneling between 
surface modes localized on the two different surfaces. The 
resonant condition corresponds to the case when the de- 
nominator in the integrand of (42) is small. Close to the 
resonance we can use the approximation 



where uji is a constant and ujq is determined by the equa- 
tion Re[e(a;o) -|- 1] = 0. In this case the heat current is 
determined bv[25j| 

Jo ~ [no(wo) - IIi(wo)] , 



where ^l « \[og{2ujah)f / {^tt). If we write Ti = Tq - AT 
and assume AT/Tq << 1 we get Jq = a AT with 



fcB7 ??^exp(7?) 
c?^ [exp(?7) — 1]^ 



(43) 



where rj — TiuJo/k-oTQ. 

Resonant photon tunnehng enhancement of the heat 
transfer is possible for two semiconductor or insula- 
tor surfaces which can support low-frequency surface 
phonon-polariton modes in the mid-infrared frequency 
region. As an example, consider two clean surfaces of 
(amorphous) silicon dioxide (Si02). The optical proper- 
ties of this material can be described using an oscillator 
modelfill 



The frequency dependent term in this expression is due 
to optical phonon's. The values for the parameters Coo, 
{a,0Ja,^a) and (6, Wfc,7&) are given in Ref. [sH- In Fig. 
13] we show the calculated heat current per unit area, Jq, 
as a function of the temperature difference AT. The 
solids have flat surfaces separated by d = 1 nm. One 
solid is at the temperature T = 296 K and the other at 
T + AT. When AT « T, the heat transfer depends 
(nearly) linearly on the temperature difference AT (see 
Fig. [5]), and we can define the heat transfer coefRcient 
= ^o/AT. In the present case (for c? = do = 1 nm) 
a — ao ^ 2x 10^ W/m^K. If the surfaces are not smooth 
but if roughness occur so that the separation d varies 
with the coordinate x = (a;, y) we have to first order in 
the expansion (41): 



/? = aQ{{do/dY 



(44) 



where (..) stands for ensemble average, or average over 
the whole surface area, and where ao is the heat transfer 
between flat surfaces separated by d = do. 

In the preset case the heat transfer is associated with 
thermally excited optical (surface) phonon's. That is, the 
electric field of a thermally excited optical phonon in one 
solid excites an optical phonon in the other solid, leading 
to energy transfer. The excitation transfer occur in both 
directions but if one solid is hotter than the other, there 
will be a net transfer of energy from the hotter to the 
colder solid. For metals, low-energy excited electron-hole 
pairs will also contribute to the energy transfer, but for 
good metals the screening of the fluctuating electric field 
by the conduction electrons leads to very ineffective heat 
transfer. However, if the metals are covered with metal 
oxide layers, and if the separation between the solids is 
smaller than the oxide layer thickness, the energy trans- 
fer may again be due mainly to the optical phonon's of 
the oxide, and the magnitude of the heat current will 
be similar to what we calculated above for (amorphous) 
silicon dioxide. 



Let us consider a high-tech application. Consider a 
MEMS device involving very smooth (amorphous) silicon 
dioxide slabs. Consider, for example, a very thin silicon 
dioxide slab rotating on a silicon dioxide substrate. Dur- 
ing operation a large amount of frictional energy may 
be generated at the interface. Assume that the disk is 
pressed against the substrate with the nominal stress or 
pressure po- This does not need to be an external applied 
force but may be due to the long-ranged van der Waals 
attraction between the solids, or due to capillary bridges 
formed in the vicinity of the (asperity) contact regions 
between the solids. The heat transfer due to the area 
of real contact (assuming purely elastic deformation) can 
be calculated from (25). Let us make a very rough esti- 
mate: Surfaces used in MEMS application have typically 
a roughness of order a few nanometers. Thus, ito ~ 1 nm 
and for (amorphous) silicon dioxide the heat conductivity 
K « 1 W/Km. Thus from (32): 



{po/E) X 10^ W/m^K 



(45) 



In a typical case the nominal pressure po may be (due to 
the van der Waals interaction and capillary bridges) be- 
tween 10^-10'^ Pa and with E « 10" Pa we get from (45) 
a « 10^ — 10^ W/Km^. If the root-mean-square rough- 
ness is of order ~ 1 nm we expect the average separation 
between the surfaces to be of order a few nanometer so 
that ((do/d)^) « 0.1 giving the non-contact contribution 
to a from the electromagnetic field of order [from (44)] 
10^ W/Km^, i.e., larger than or of similar magnitude as 
the contribution from the area of real contact. 

(b) Contribution to a from heat transfer via the 
surrounding gas or liquid 

Consider two solids with flat surfaces separated by a 
distance d. Assume that the solids are surrounded by a 
gas. Let A be the gas mean free path. If d >> A the 
heat transfer between the solids occurs via heat diffusion 
in the gas. If d << A the heat transfer occurs by bal- 
listic propagation of gas molecules from one surface to 
the other. In this case gas molecules reflected from the 
hotter surface will have (on the average) higher kinetic 
energy that the gas molecules reflected from the colder 
surface. This will result in heat transfer from the hotter 
to the colder surface. The heat current is approximately 
given bv[3^ 



Jo 



,AT 



d+ aA 



where a is a number of order unity and which depend on 
the interaction between the gas molecules and the solid 
walls For air (and most other gases) at the normal at- 
mospheric pressure and at room temperature A « 65 nm 
and Kgas ~ 0.02 W/mK. For contacting surfaces with 
surface roughness we get to first order in the expansion 
in (41): 



Kgas((d-HA) 1) 



(46) 



where (..) stand for ensemble average or averaging over 
the surface area. Eq. (46) also holds if the surfaces are 
surrounded by a liquid rather than a gas. In this case 
Kgas must be replaced with the liquid heat conductivity 
Kiiq and in most cases one can put A equal to zero. 

If we again consider a MEMS application where the 
average surface separation is of order nm we can neglect 
the d-dependence in (46) and get a ~ Kgas/A « 3 x 
10^ W/m^K which is similar to the contribution from 
the electromagnetic coupling. 

(c) Contribution to a from heat transfer via 
capillary bridges 

If the solid walls are wet by water, in a humid at- 
mosphere capillary bridges will form spontaneous at the 
interface in the vicinity of the asperity contact regions. 
For very smooth surfaces, such as in MEMS applications, 
the fluid (in this case water) may occupy a large region 
between the surfaces and will then dominate the heat 
transfer between the solids. Similarly, contamination lay- 
ers (mainly organic molecules) which cover most natural 
surfaces may form capillary bridges between the contact- 
ing solids, and contribute in an important way to the 
heat transfer coefficient. The fraction of the interfacial 
surface area occupied by fluid bridges, and the separa- 
tion between the solids in the fluid covered region, can 
be calculated using the theory developed in Ref. [11]. 
From this one can calculate the contribution to the heat 
transfer using (46): 
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where P(u) is the distribution of interfacial separation u, 
and ^0 the nominal contact area. The lower cut-off a in 
the integral is a distance of order a molecular length and 
dx is the maximum height of the liquid bridge which, 
for a system in thermal equilibrium and for a wetting 
liquid, is of order the Kelvin length. Note that P(u) is 
normalized and that 



dK 



duAoP{u) = AA 



(48) 



is the surface area (projected on the xy-plane) where the 
surface separation is between a < m < c?k- 

3. Contact mechanics: short review and basic 
equations 

The theory of heat transfer presented above depends 
on quantities which can be calculated using contact me- 
chanics theories. Thus, the heat flux through the non- 
contact area (Sec. 2.2.2) depends on the average of some 
function f[d{x.)] of the interfacial separation d{x.). If 
P{u) denote the probability distribution of interfacial 
separation u then 

/•OO 

{fid)) = / du f{u)P{u) (49) 

J a 




FIG. 5: An rubber block (dotted area) in adhesive contact 
with a hard rough substrate (dashed area). The substrate 
has roughness on many different length scales and the rubber 
makes partial contact with the substrate on all length scales. 
When a contact area is studied at low magnification it appears 
as if complete contact occur, but when the magnification is 
increased it is observed that in reality only partial contact 



where a is a short-distance cut-off (typically of molecu- 
lar dimension). The contribution from the area of real 
contact depends on the elastic energy Uc\ stored in the 
asperity contact regions [see Eq. (23)]. In the limit of 
small contact pressure Uc\ = PqUq, where uq is a length 
which is of order the root-mean-square roughness of the 
combined roughness profile. All the quantities P{u), Uei 
and uq can be calculated with good accuracy using the 
contact mechanics model of Persson. Here we will briefly 
review this theory and give the basic equations relevant 
for heat transfer. 

Consider the frictionless contact between two elastic 
solids with the Young's elastic modulus Eq and Ei and 
the Poisson ratios vq and vi . Assume that the solid sur- 
faces have the height profiles /io(x) and /ii(x), respec- 
tively. The elastic contact mechanics for the solids is 
equivalent to those of a rigid substrate with the height 
profile /i(x) = /io(x) -I- /ii(x) and a second elastic solid 
with a flat surface and with the Young's modulus E and 
the Poisson ratio i' chosen so that 34] 
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(50) 



The contact mechanics formalism developed 
elsewhere [13, U, [12] is based on the studying 
the interface between two contacting solids at differ- 
ent magnification C. When the system is studied at 
the magnification C it appears as if the contact area 
(projected on the xy-plane) equals A(C,), but when the 
magnification increases it is observed that the contact 
is incomplete, and the surfaces in the apparent contact 
area A{C,) are in fact separated by the average distance 
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FIG. 6: An asperity contact region observed at the magnifica- 
tion It appears that complete contact occur in the asperity 
contact region, but when the magnification is increasing to the 
highest (atomic scale) magnification , it is observed that the 
solids are actually separated by the average distance ii(C)- 



m(C), see Fig. [51 The (apparent) relative contact area 
A{Q/Aq at the magnification ( is given byjol. ITU 
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where 

G(C) = - 
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where the surface roughness power spectrum 
1 



C{q) 



d^x {h{-K)h{Q))e- 



(51) 
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where (...) stands for ensemble average. The height pro- 
file /i(x) of the rough surface can be measured routinely 
today on all relevant length scales using optical and sty- 
lus experiments. 

We define ui(C) to be the (average) height separating 
the surfaces which appear to come into contact when the 
magnification decreases from C to C — AC, where AC is a 
small (infinitesimal) change in the magnification. ui(C) 
is a monotonically decreasing function of C, and can be 
calculated from the average interfacial separation m(C) 
and A(C) using (see Ref. (ll| ) 



where 11 1 



^.i(C) = s(C) + "'(C)^(C)M'(C), 



u{0 = / dq q^C{q)w{q) 



(54) 
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FIG. 7: An elastic block squeezed against a rigid rough sub- 
strate. The separation between the average plane of the sub- 
strate and the average plane of the lower surface of the block 
is denoted by u. Elastic energy is stored in the block in the 
vicinity of the asperity contact regions. 



where E* = E/{1 - v"^), and where p{C) — po^oM(C) 
and 

r /2 

MqX) = dq' q''C{q')^ 

The distribution of interfacial separations 

P{u) = (S[u-u{x)]) 

where ■(i(x) = (i(x) is the separation between the surfaces 
at point X. As shown in Ref. [llj we have (approxi- 
mately) 

/oc 
dc [-A'iCMu-u.iO] 



Thus we can write (49) as 



(/(d)) = / dC [-A'(C)]/[^i(C)] 



(56) 



where Ci is defined by wi(Ci) = o- 

Finally, the elastic energy Ud and the length parameter 
Mo can be calculated as follows. The elastic energy Uc\ 
has been studied in Ref. Ml: 



U.i^AoE*- 



dq q^W{q,p)C{q). 



(57) 



In the simplest case one take W{q,p) = P{q,p) = 
A{()/Aq is the relative contact area when the interface is 
studied at the magnification C = q/qo, which depends on 
the applied pressure p — po- A more accurate expression 
is 



W{q,p) = P{q,p) [7 + (1 - l)P^{q,p)] 



(58) 
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FIG. 8: Surface roughness power spectrum C{q) as a function 
of the wavevector g on a log-log scale (with 10 as basis) . For 
a typical surface used in MEMS applications with the root 
mean square roughness 2.5 nm when measured over an area 
10 nm X 10 nm. 



However, in this case one also need to modify (55) appro- 
priately (see Ref. [nj). The parameter 7 in (58) seams 
to depend on the surface roughness. For self-afRne frac- 
tal surfaces with the fractal dimension Df « 2.2 we have 
found that 7 « 0.5 gives good agreement between the 
theory and numerical studies As Df ^ 2 analysis of 
numerical data indicate that 7 — > 1. 



For small pressures one can show that[lll|: 



p = f3E*e- 
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where 



uo = V^^ / d,q q C{q)'w{q), 
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where w(g) = w{q, 1), and where 
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where (for 7 = 1) e = 0.7493. 
4. Numerical results 

In this section we present numerical results to illus- 
trate the theory. We focus on a MEMS-like application. 
In Fig. [8] we show the surface roughness power spec- 
trum C{q) as a function of the wavevector g on a log- 
log scale (with 10 as basis) for a typical surface used in 
MEMS applications, with the root mean square rough- 
ness 2.5 nm when measured over an area 10 /im x 10 /im. 
In Fig. [9] we show for this case the contribution to the 
heat transfer coefficient a from the direct contact area, 
and the non-contact contribution due to the fluctuat- 
ing electromagnetic (EM) field and due to heat transfer 
via the surrounding gas. In the calculation of the EM- 
contribution we have used (44) with ag = 2.0 MW/m^K 
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FIG. 9: The contribution to the heat transfer coefficient a 
from the direct contact area, and the non-contact contribution 
due to the fluctuating electromagnetic (EM) field and due 
to heat transfer via the surrounding gas. For a randomly 
rough surface with the (combined) surface roughness power 
spectrum shown in Fig. [S] 
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FIG. 10: The logarithm (with 10 as basis) of the contribution 
to the heat transfer coefficient a from the real contact areas, 
and from the water in the capillary bridg function 
of the relative (water) humidity. For a randomly rough sur- 
face with the (combined) surface roughness power spectrum 
shown in Fig. [S] The squeezing pressure po = 4 MPa and 
the effective solid elastic modulus E* — 86 GPa. The heat 
conductivity of water Kfluid = 0.58 W/mK. 



(and do = 1 nm). For the contribution from the sur- 
rounding gas we have used (46) with Kgas = 0.024 W/mK 
and A = 65 nm (and a ~ 1). For the contact contribu- 
tion we used (25) with k — 1 W/mK. In all calculations 
we have assumed E* — 86 GPa and that the contact is 
elastic (no plastic yielding). 

We have also studied the contribution to the heat 
transfer from capillary bridges which on hydrophilic sur- 
faces form spontaneous in a humid atmosphere. The cap- 
illary bridges gives an attractive force (to be added to the 
external squeezing force) , which pulls the solids closer to- 



gether. We have used the theory presented in Ref. [33| to 
include the influence of capillary bridges on the contact 
mechanics, and to determine the fraction of the interface 
area filled with fluid at any given relative humidity. In 
Fig. [TO] we show the logarithm (with 10 as basis) of the 
contribution to the heat transfer coefficient a from the 
real contact areas, and from the water in the capillary 
bridg function of the relative (water) humidity. 

For relative humidity below ~ 0.4 the contribution to 
the heat transfer from capillary bridges decreases roughly 
linearly with decreasing humidity (and vanish at zero hu- 
midity), and for relative humidity below ^ 0.015 the heat 
transfer via the area of real contact will be more impor- 
tant than the contribution from the capillary bridges. 
However the contribution from heat transfer via the air 
or vapor phase (not shown) is about ^ 0.3 MW/m^K 
(see Fig. |9]), and will hence give the dominant contribu- 
tion to the heat transfer for relative humidity below 0.3. 
The small increase in the contribution from the area of 
real contact for relative humidity around ^ 0.94 is due 
to the increase in the contact area due to the force from 
the capillary bridges. For soft elastic solids (such as rub- 
ber) this effect is much more important: see Ref. [33| for 
a detailed discussion of this effect, which will also affect 
(increase) the heat transfer in a drastic way. 

We note that heat transfer via capillary bridges 
has recently been observed in nanoscale point contact 
experiments |35j. In this study the authors investigated 
the heat transfer mechanisms at a ~ 100 nm diameter 
point contact between a sample and a probe tip of a 
scanning thermal microscope. They observed heat trans- 
fer both due to the surrounding (atmospheric) air, and 
via capillary bridges. 

5. Role of adhesion and plastic deformation 

In the theory above we have assumed that the solids 
deform purely elastically. However, in many practical sit- 
uations the solids will deform plastically at short enough 
length scale. Similarly, in many practical situations, in 
particular for elastically soft solids, the area of real con- 
tact may depend strongly on the adhesive interaction 
across the contacting interface. Here we will briefly dis- 
cuss under which circumstances this will affect the heat 
transfer between the solids. 

The contribution to the heat transfer from the area 
of real contact between two solids depends on the elas- 
tic energy Ue\ stored in the asperity contact regions, or, 
at small enough applied loads, on the length parame- 
ter uq. For most randomly rough surfaces these quanti- 
ties are determined mainly by the long- wavelength, large 
amplitude surface roughness components. Similarly, the 
interfacial separation, which determines the non-contact 
contribution to the heat transfer, depends mainly on the 
long-wavelength, large amplitude surface roughness com- 
ponents. On the other hand, plastic deformation and ad- 
hesion often manifest themself only at short length scales. 
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FIG. 11: The elastic Aei and plastic contact area as a func- 
tion of magnification on a log- log scale (with 10 as basis). The 
penetration hardness ~ 4 GPa and the applied pressure 
po — 4 MPa. Also shown is the asperity-induced elastic en- 
ergy (/ei(C) in units of the full elastic energy C/ei(Ci) obtained 
when all the roughness (with wavevectors below gi = Ci9o) is 
included. The vertical dashed line indicate the magnification 
where A^i = Api. 



corresponding to high magniflcation. For this reason, 
in many cases one may assume purely elastic deforma- 
tion when calculating the heat transfer, even if, at short 
enough length scale, all asperities have yielded plasti- 
cally, or the adhesion has strongly increased the (appar- 
ent) contact area. Let us illustrate this with the amor- 
phous silicon dioxide system studied in Sec. 4. 

In Fig. [11] we show the elastic and plastic contact area 
as a function of magnification on a log-log scale (with 
10 as basis). Also shown is the asperity-induced elastic 
energy Uei{C) in units of the full elastic energy C/ei(Ci) ob- 
tained when all the roughness (with wavevectors below 
qi — Ciqo) is included. Note that about 90% of the full 
elastic energy is already obtained at the magnification 
where the elastic and plastic contact areas are equal, and 
about 60% of the full elastic energy is obtained when 
^piMci ~ 0.01. Thus, in the present case, to a good 
approximation, we can neglect the plastic deformation 
when studying the heat transfer. In the calculation we 
have assumed the penetration hardness ay = 4 GPa and 
the squeezing pressure po = 4 MPa. Thus, at high magni- 
fication, where all the contact regions are plastically de- 
formed, the relative contact area A/Aq — po/ay = 0.001 
in good agreement with the numerical data in Fig. [TT] 

If necessary, it is easy to include adhesion and plastic 
deformation when calculating the heat transfer coefficient 
a. Thus (26b) is also valid when adhesion is included, at 
least as long as adhesion is treated as a contact interac- 
tion. However, in this case the interfacial stiffness dpp/du 
must be calculated including the adhesion (see Ref. [36(). 
Plastic deformation can be included in an approximate 
way as follows. If two solids are squeezed together at 




FIG. 12: Temperature distribution of rubber tread (thickness 
d) in contact with the air. The air temperature (for z > d) 
and the temperature at the outer {z = d) and inner {z — 0) 
rubber surfaces are denoted by Tair, Ti and To, respectively. 



the pressure po they will deform elastically and, at short 
enough length scale, plastically. If the contact is now re- 
moved the surfaces will be locally plastically deformed. 
Assume now that the surfaces are moved into contact 
again at exactly the same position as the original con- 
tact, and with the same squeezing pressure po applied. 
In this case the solids will deform purely elastically and 
the theory outlined in this paper can be (approximately) 
applied assuming that the surface roughness power spec- 
trum C{q) of the (plastically) deformed surface is known. 
In Ref. 12] we have described an approximately way of 
how to obtain C{q) from C{q) by defining (with q = (qo) 
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where = Fn/cty- The basic picture behind this def- 
inition is that surface roughness at short length scales 
get smoothed out by plastic deformation, resulting in an 
effective cut-off of the power spectrum for large wavevec- 
tors (corresponding to short distances). 

6. Application to tires 

Here we will briefly discuss heat transfer in the context 
of tires. The rolling resistance /^r of a tire determines the 
heat production in a tire during driving on a strait planar 
road at a constant velocity v. In a stationary state the 
energy produced per unit time, W = /j,r,Fnw, must equal 
the transfer of energy per unit time, from the tire to the 
surrounding atmosphere and to the road surface. Here 
we will briefly discuss the relative importance of these 
two different contributions to the heat transfer. 



Assume for simplicity that the frictional heat is pro- 
duced uniformly in the tread rubber, and assume a tire 
without tread pattern. Let z be a coordinate axis perpen- 
dicular to the rubber surface. In this case at stationary 
condition the temperature in the tread rubber satisfies 
T"{z) = —q/n where q is the frictional heat produced 
per unit volume and unit time. We assume that the heat 
current vanish at the inner rubber surface (z — 0, see Fig. 
ini), so that r'(0) = 0. Thus we get T[z) = Tq - qz^/2K. 
The heat current at the outer rubber surface 

Jo = -KT'{d) = qd. (62) 

The temperature of the outer surface of the tread rubber 

Ti = T{d) = To- qd^/2K (63) 

Let us now assume that the heat transfer to the surround- 
ing 



Jo ^ a{Ti- Tair) 

Combining (62)-(64) gives 



Ti = To - 



Jb — Tair 

1 + 2K/da 



(64) 



(65) 



For rubber k « 0.2 W/mK and with d — 1 cm and 
a « 100 W/m^K, as is typical for (forced) convective 
heat transfer between a tire and (dry) air (see Appendix 
E and Ref. [stI), we get 

Ti « 0.3To + O.TTair. 

The temperature profile is shown (schematically) in Fig. 
[T^ In reality, the heat production, even during pure 
rolling, will be somewhat larger close to the outer sur- 
face of the tread and the resulting temperature profile 
in the tread rubber will therefore be more uniform than 
indicated by the analysis above. 

Let us now discuss the relative importance of the con- 
tributions to the heat transfer to the air and to the road. 
We assume that the heat transfer to the atmosphere and 
to the road are proportional to the temperature differ- 
ence Ti — Tair and Ti — Tj-oad , respectively. We get 

AiR-pNt' = aait-4surf (^i - Tair) -|- aroad^O (Tl - Tmad) (66) 

where ^surf is the outer surface area of the tread, and 
Aq the nominal tire-road footprint area. For rubber in 
contact with a road surface k in Eq. (22) is ~ 0.2 W/mK 
and with po/E* « 0.04 and uq ~ 10"'^ m (as calculated 
for a typical case) we get aroad ~ 10 W/m^K which is 
smaller than the contribution from the forced convection. 
Since the nominal contact area between the tire and the 
road is much smaller than the total rubber tread area, 
we conclude that the contribution from the area of real 
contact between the road and the tire is rather unimpor- 
tant. During fast acceleration wear process may occur, 



involving the transfer of hot rubber particles to the road 
surface, but such processes will not considered here. In 
addition, at the inlet of the tire-road footprint area, air 
may be be compressed and then rapidly squeezed out 
from the tire-road contact area resulting in strong forced 
convective cooling of the rubber surface in the contact 
area. A similar process involving the inflow of air occur 
at the exit of the tire-road footprint area. A detailed 
study of this complex process is necessary in order to ac- 
curately determine the heat transfer from a tire to the 
surrounding atmosphere and the road surface. 

For a passenger car tire during driving on a strait pla- 
nar road at a constant velocity u, the tire temperature 
which follows from (66) is in reasonably agreement with 
experiment. Thus, using (66) we get 



heating mixing 



AT = Ti - Tair 



(67) 



and with aair — 100 W/ni^K, Agm-f ~ 0.5 vc? and « 
0.02, Fn = 3500 N and w = 30 m/s we get AT k. 40 °C. 

The discussion above has focused on the stationary 
state where the heat energy produced in the tire per unit 
time is equal to the energy given off to the surrounding 
per unit time. However, for a rolling tire it may take a 
very long time to arrive at this stationary state. In the 
simplest picture, assuming a uniform temperature in the 
tire rubber, we get from energy conservation 

=^-^(T-Ta.) 

or, if r(0) = Tair, 

a V 



where the relaxation time r — pC\d/a ~ 200 sec. In 
reality, the temperature in the tire is not uniform, and 
this will introduce another relaxation time r', defined as 
the time it takes for heat to diffuse a distance d, which 
is of order r' = pCyd"^ /k. The ratio t'/t = ad/n. For 
rubber k « 0.2 W/mK and assuming d — 1 cm and a — 
100 W/m^K gives t'/t sa 5 or r' sa 10^ sec. Experiment 
have shown that it typically takes ~ 30 minutes to fully 
build up the tire temperature during rolling [STj. 

Rubber friction depends sensitively on the tempera- 
ture of the rubber, in particular the temperature close to 
the rubber surface in contact with the road. The tem- 
perature in the surface region of a tire varies rapidly in 
space and time, which must be considered when calculat- 
ing the rubber friction[lj. The shortest time and length 
scales are related to the contact between the road as- 
perities and the rubber surface in the tire-road footprint 
contact area. During slip this generate intense heating 
which varies over length scales from a few micrometer to 
several mm, and over time scales shorter than the time a 
rubber patch stays in the footprint, which typically may 
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FIG. 13: Experiment to test the theory predictions for the 
heat transfer across interfaces. The increase in the tempera- 
ture Ti it) of the water in the lower container, with increasing 
time t, determines the heat transfer between the upper and 
lower water container. 



be of order a few milliseconds. During this short time 
very little heat is transferred to the surrounding, and 
very little heat conduction has occurs inside the rubber, 
i.e., the heat energy mainly stays where it is produced 
by the internal friction in the rubber. This result in a 
flash temperature effect, which has a crucial influence on 
rubber friction[l|. However, rubber friction also depends 
on the background temperature (usually denoted by Tg), 
which varies relatively slowly in space and time, e.g., on 
time scales from the time ~ 0.1 s it takes for the tire to 
perform a few rotations, up to the time ~ 30 minutes 
necessary to build up the full tire temperature after any 
change in the driving condition (e.g., from the start of 
driving). Note that the time variation of the background 
temperature Tq depends on the surrounding (e.g., the air 
and road temperatures, humidity, rain, ...) and on the 
driving history, while the flash temperature effect mainly 
depends on the slip history of a tread block (or rubber 
surface patch) in the footprint contact area, but not on 
the outside air or road temperature, or atmospheric con- 
dition. 

7. A new experiment 

We have performed a very simple experiment to test 
the theoretical predictions for the heat transfer. The 
setup consists of two containers, both filled with dis- 
tilled water, standing on top of each other with a thin 
silicon rubber film in between. The upper container is 
made from copper (inner diameter 5 cm), and the water 



is heated to the boiling temperature (i.e., Tq = 100°C). 
The lower container is made from PMMA with a cylin- 
drical copper block at the top. To study the effect of 
surface roughness on the heat transfer, the copper block 
can be replaced by another copper block with different 
surface roughness. In the experiments presented below 
we used 3 copper blocks with different surface roughness. 

The temperature Ti (t) of the water in the lower con- 
tainer will increase with time t due to the heat current Jo 
flowing from the upper container to the lower container: 



Jo = pCyTid 



(68) 



where d is the height of the water column in the lower 
container (in our experiment d = 3.5 cm), and where p 
and Cy are the water mass density and heat capacity 
respectively. We measure the temperature of the water 
in the lower container as a function of time, starting at 
25°C. To obtain a uniform temperature of the water in 
the lower container we mix it using a (magnetic-driven) 
rotating metal bar. 

We have investigated the heat transfer using copper 
blocks with different surface roughness. To prepare the 
rough surfaces, we have pressed annealed (plastically 
soft) copper blocks with smooth surface against sandpa- 
per, using a hydraulic press. We repeated this procedure 
several times to obtain randomly rough surfaces. The 
roughness of the copper surfaces can be changed by us- 
ing sandpaper of different grade (consisting of particles 
with diff'erent (average) diameter). Due to the surface 
roughness, the contact between the top surface of the 
lower container and the thin silicon rubber sheet (thick- 
ness do = 2.5 mm) attached to the upper container, is 
only partial. The bottom surface of the upper container 
has been highly polished and we can neglect the heat re- 
sistance at this rubber-copper interface. Thus, most of 
the resistance to the heat flow arises from the heat dif- 
fusion through the rubber sheet, and from the resistance 
to the heat flow at the interface between the rubber and 
the rough copper block. 

The rubber sheet (elastic modulus E = 2.5 MPa, 
Poisson ration u = 0.5) was made from a silicone elas- 
tomer (PDMS). We have used Polydimcthylsiloxane be- 
cause of its almost purely elastic behavior on the time 
scales involved in our experiments. The PDMS sam- 
ple was prepared using a two-component kit (Sylgard 
184) purchased from Dow Corning (Midland, MI). This 
kit consists of a base (vinyl-terminated polydimethyl- 
siloxane) and a curing agent (mcthylhydrosiloxane- 
dimethylsiloxane copolymer) with a suitable catalyst. 
Prom these two components we prepared a mixture of 
10:1 (base/cross linker) in weight. The mixture was de- 
gassed to remove the trapped air induced by stirring from 
the mixing process and then poured into cylindrical casts 
(diameter 5 cm and height = 2.5 mm). The bottom of 
these casts were made from glass to obtain smooth sur- 



faces (negligible roughness). The samples were cured in 
an oven at 80 °C for over 12 hours. 
Using (3) we can write 



Jo 



To-T,{t) 
doKQ^ + a-^ 



(69) 



where no the heat conductivity of the rubber. Here 
we have neglected the influence of the copper blocks on 
the heat transfer resistance, which is a good approxima- 
tion because of the high thermal conductivity of copper. 
Combining (68) and (69) gives 

Tfi=To-Ti{t) 



where the relaxation time 



To = pCyd 



do , 1 



If we assume that to is time independent, we get 

Ti (i) = To + [Ti (0) - To] e-*/^" . (70a) 

In the study above we have assumed that there is no 

heat transfer from the lower container to the surrounding. 
However, if necessary one can easily take into account 
such a heat transfer: If we assume that the heat transfer 
depends linearly on the temperature difference between 
the water and the surrounding we can write 

Jl = ai (Ti - Tsurr) 

In this case it is easy to show that (70a) is replaced with 

Ti(i) =Ta + [Ti(0) - Ta] e-*/^ (705) 

where Ta is the temperature in the water after a long 
time (stationary state where Jo = Ji), and where the 
relaxation time t now is given by 



T = pCyd 



Ta Tguj-j- / do 1 



To — Tsurr V '^0 CC 



The heat transfer across the rubber-copper interface 
can occur via the area of real contact, or via the non- 
contact area via heat diffusion in the thin air film or 
via radiative heat transfer. Since all these heat transfer 
processes act in parallel we have 



a 



Let us estimate the relative importance of these dif- 
ferent contributions to a. Using the (diffusive) heat 
conductivity of air Kgas « 0.02 W/mK and assuming 
(d-^) = (20 /im)-i gives 

agas = Kgas((d + A)-i) « Kg^{d-^) « 1000 W/m^K. 
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FIG. 14: The surface roughness power spectrum of the three 
copper surfaces used in the experiment. The surfaces 1, 2 and 
3 have the root-mean-square roughness 42, 88 and 114 /im, 
respectively. 



Let US assume that po « 0.01 MPa, E* ^ 2 MPa, uq 
10 and (for rubber) kq = 0.2 W/mK. Thus 



E*uo 



100 W/m^K. 



Here we have used that k w kq (since the heat conduc- 
tivity Ki of copper is much higher than for the rubber). 
FinaUy, assuming the radiative heat transfer is weU ap- 
proximated by the Stefan-Boltzmann law and assuming 
that (To - Ti)/Ti « 1, we get with Tq = 373 K 
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FIG. 15: The variation of the of the heat transfer coefficient 
from the contact area (aeon) and from the air-gap (cvgas) with 
the squeezing pressure. The surfaces 1, 2 and 3 have the 
power spectra's shown in Fig. 1141 



ttrad 



10 W/m^K 



Note that arad is independent of the squeezing pressure 
po, while aeon ~ Po- The pressure dependence of agas will 
be discussed below. 

In the experiment reported on below the silicon rubber 
film has the thickness do 
100 W/m^K. Thus 



2.5 mm so that dg ^kq 



1 



1 

a 



1 

Too 



1 



1000+ 100+ 10 



(W/m^K) 



and it is clear from this equation that in the present case 
the thin rubber film will give the dominant contribution 
to the heat resistance. This is in accordance with our 
experimental data presented below. 

8. Experimental results and discussion 

To test the theory we have performed the experiment 
described in Sec. 7. We have performed experiments 
on four different (copper) substrate surfaces, namely one 
highly polished surface (surface 0) with the root-mean- 
square (rms) roughness 64 nm, and for three rough sur- 
faces with the rms roughness 42, 88 and 114 /im. In Fig. 
[T3]we show the surface roughness power spectrum of the 
three latter surfaces. Including only the roughness with 



wavelength above ~ 30 /um, the rms slope of all three 
surfaces are of order unity, and the normalized surface 
area A/Aq « 1.5 in all cases. 

In Fig. [15] we show for the surfaces 1, 2 and 3, the 
pressure dependence of heat transfer coefficient from the 
contact area (aeon) and from the air-gap (agas). Note 
that both aeon and agas varies (nearly) linearly with po- 
The latter may at first appear remarkable because we 
know that at the low (nominal) squeezing pressures used 
in the present calculation (where the area of real contact 
varies linearly with pq), the average surface separation 
u = (u) depends logarithmically on po- However, the 
heat transfer via heat diffusion in the air gap depends 
on ((u + A)^^) which depends on po almost linearly as 
long as M >> A, which is obeyed in our case. This can 
be understood as follows: (u) is determined mainly by 
the surface regions where the surface separation is close 
to its largest value. On the other hand {{u + A)~^) is 
determined mainly by the surface regions where u is very 
small, i.e., narrow strips (which we will refer to as bound- 
ary strips) of surface area close to the area of real con- 
tact. Now, for small po the area of real contact increases 
linearly with po while the distribution of sizes of the con- 
tact regions is independent of pq. It follows that the 
total area of the boundary strips will also increase lin- 




separation (mm) 

FIG. 16: The variation of the cumulative probability with 
the height (or gap-separation) u. The surfaces 1 and 3 (top) 
and 2 (bottom) have the power spectra's shown in Fig. 1141 
For each surface the curves are for the nominal squeezing 
pressures (from left to right): 11.8, 23.7, 35.5, 47.3, 59.2 and 
71.0 kPa. 

early with pq. Thus, since {{u + A)^^) is determined 
mainly by this surface area, it follows that {{u + A)^^) 
will be nearly proportional to pq. We note that in Fig. [5] 
Qfgas is nearly pressure independent, but this is due to the 
fact that the (combined) surface in this case is extremely 
smooth (root- mean-square roughness 2.5 nm) so that the 
M-terni in {(u + A)~^) can be neglected compared to the 
gas mean free path A, giving a nearly pressure indepen- 
dent gas heat transfer coefRcient. However, in the system 
studied above u is much larger than A and the result is 
nearly independent of A. 

Note that in the present case (see Fig. [15]) Ofgas >> 
Q!con so that the present experiment mainly test the the- 
ory for the heat flow in the air gap. 

In Fig. [16] we show the variation of the cumulative 
probability with the height (or gap-separation) u for the 
surfaces 1 and 3 (top) and 3 (bottom). 

In Fig. [17] we show the measured (dots) and cal- 
culated [using (70b)] (solid lines) temperature in the 
lower container as a function of time. Results are for 
all four surfaces and for the nominal squeezing pressure 
Po = 0.012 MPa. In Fig. [TH]we show the measured (dots) 



and calculated (solid lines) temperature in the lower con- 
tainer as a function of time. Results are for surface 2 for 
the nominal squeezing pressure po = 0.012 (lower curve) 
and 0.071 MPa (upper curve). Note that there is no fit- 
ting parameter in the theory calculations, and the agree- 
ment between theory and experiment is relative good. 

The heat resistance of the system studied above is dom- 
inated by the thin rubber film. The reason for this is the 
low heat conductivity of rubber (roughly 100 times lower 
than for metals). For direct metal-metal contact the con- 
tact resistance will be much more important. However, 
for very rough surfaces it is likely that plastic flow is ob- 
served already at such low magniflcation (corresponding 
to large length scales) that it will affect the contact re- 
sistance. Nevertheless, it is interesting to compare the 
theory predictions for elastic contact with experimental 
data for metal-metal contacts. 

In Fig. [19] we show the measured heat transfer coef- 
flcient for metal-metal contacts with steel, copper and 
aluminumlil]. The surfaces have the effective (or com- 
bined) rms surface roughness /ij-ms = 7.2 pm (steel), 
2.2 fiui (Cu) and 5.0 pm (Al). Assume that the variation 
of a with Po is mainly due to the area of real contact, 
i.e., we neglect the heat transfer via the thin air fllm be- 
tween the surfaces. Fitting the data points in Fig. [19] 
with strait lines gives the slope da/c?po(exp) (in units of 
m/sK): 

2 X lO""* (steel), 7 x 10"^ (Cu), 1.2 x 10"^ (Al) 

Using (26a) with uq w OAhj-^s (here we have assumed 
7 = 0.4) gives da/dpo (theory) = k/E*uo: 

1 X lO""* (steel), 4 x 10"^ (Cu), 1.3 x 10"^ (Al) 

The agreement between theory and experiment is very 
good taking into account that plastic deformation may 
have some influence on the result, and that an accurate 
analysis requires the full surface roughness power spec- 
trum C (q) (in order to calculate uq accurately, and in or- 
der to include plastic deformation if necessary (see Sec. 
5)), which was not reported on in Ref. [s^. We note 
that experimental results such as those presented in Fig. 
1191 are usually analyzed with a phenomenological model 
which assumes plastic flow and neglect elastic deforma- 
tion. In this theory the heat transfer coefficient [s^ 



is proportional to the rms surface slope s, but it is well 
known that this quantity is dominated by the very short- 
est wavelength roughness which in fact makes the theory 
ill-deflned. In Ref. [11] the data presented in Fig. [T51 
was analyzed using (71) with s = 0.035, 0.006 and 0.03 
for the steel, Cu and Au surfaces, respectively. However, 
analysis of polished surfaces with similar rms roughness 
as used in the experiments usually gives slopes of order 
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FIG. 17: The measured (dots) and calculated (solid lines) 
temperature in the lower container as a function of time. Re- 
sults are for all four surfaces and for the nominal squeezing 
pressure po = 0.012 MPa. 
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FIG. 18: The measured (dots) and calculated (solid lines) 
temperature in the lower container as a function of time. 
Results are for surface 2 for the nominal squeezing pressure 
Po = 0.012 (lower curve) and 0.071 MPa (upper curve). 



unity when all roughness down to the nanometer is in- 
cluded in the analysisfi^. Using s w 1 in (71) gives heat 
transfer coefBcients roughly 100 times larger than ob- 
served in the experiments. (In our theory [Eq. (26a)] 
s/(7y in (71) is replaced by 1/E*, and since typically 
E* I<jy ~ 100, our theory is consistent with experimen- 
tal observations.) 4l| We conclude that the theory behind 
(71) is incorrect or incomplete. A theory which includes 
both elastic and plastic deformation was described in Sec. 
5. 



9. Electric contact resistance 

It is easy to show that the problem of the electrical con- 
tact resistance is mathematically equivalent to the prob- 
lem of the thermal contact resistance. Thus, the electric 
current (per unit nominal contact area) Jq through an 
interface between solids with randomly rough surfaces 
can be related to the electric potential drop at the 
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FIG. 19: Variation of the heat transfer coefficient q with 
the squeezing pressure po for metal-metal contact with steel, 
copper and aluminum. The surfaces have the effective (or 
combined) root-mean-square surface roughness values /irms = 
7.2 /xm (steel), 2.2 ^m (copper) and 5.0 /xm (aluminum). The 
heat conductivity of the metals are k = 54 W/mK (steel), 
381 W/mK (copper) and 174 W/mK (aluminum). Based on 
experimental data from Ref. 'ss']. 



interface via Jq — a'A^ where, in analogy with (25), 



a 



(72) 



where k' is the electrical conductivity. However, from a 
practical point of view the problem of the electrical con- 
tact resistance is more complex than for the heat contact 
resistance because of the great sensitivity of the electric 
conductivity on the type of material (see Appendix D). 
Thus, in a metal-metal contact the contact resistance will 
depend sensitively on if the thin insulating oxide layers, 
which covers most metals, are fractured, so that direct 
metal-metal contact can occur. On the other hand, in 
most cases there will be a negligible contribution to the 
electric conductivity from the non-contact regions. 

10. Summary and conclusion 

We have studied the heat transfer between elastic 
solids with randomly rough but nominally flat surfaces 
squeezed in contact with the pressure po • Our approach 
is based on studying the heat flow and contact mechanics 
in wavevector space rather than real space which has the 
advantage that we do not need to consider the very com- 
plex fractal-like shape of the contact regions in real space. 
We have included both the heat flow in the area of real 
contact as well as the heat flow across the non-contact 
surface region. For the latter contribution we have in- 
cluded the heat transfer both from the fluctuating elec- 
tromagnetic field (which surrounds all material objects), 
and the heat flow via the surrounding gas or liquid. We 
have also studied the contribution to the heat transfer 
from capillary bridges, which form spontaneously in a 
humid atmosphere (e.g., as a result of organic and water 



contamination films which occur on most solid surfaces 
in the normal atmosphere). We have presented an illus- 
trative application relevant for MEMS applications in- 
volving very smooth amorphous silicon dioxide surfaces. 
In this case we find that all the mentioned heat transfer 
processes may be roughly of equal importance. 

Wc have briefly discussed the role of plastic deforma- 
tion and adhesion on the contact heat resistance. We 
have pointed out that even if plastic deformation and ad- 
hesion arc important at short length scale (or high magni- 
fication) they may have a negligible influence on the heat 
transfer since the elastic energy stored in the asperity 
contact regions, which mainly determines both the inter- 
facial separation and the contact heat transfer coefficient, 
is usually mainly determined by the long- wavelength sur- 
face roughness components, at least for fractal-like sur- 
faces with fractal dimension Df < 2.5 (which is typically 
obeyed for natural surfaces and surfaces of engineering 
interest). 
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Appendix A 

In Sec. 2.2.1 we have assumed that 

^ / d^gi(|AJ,(q)p) -1 / d'q-{\AaMr) (Al) 

This equation is a c;onsequenc;e of the fact that for clastic 
solids with randomly rough surfaces the heat transfer co- 
efficient depends mainly on the geometrical distribution 
of the contact area. This can be understood as follows. 
Let x„ denote the center of the contact spot n and let /„ 
be the heat current through the same contact spot. We 
now approximate 

Jz{y^) ~ ^-fn^(x-X„). 



Thus 



and 



Thus the left hand side (LHS) of (Al) becomes 



LHS 



^0 



1 iq-(xm-x„) 



(A2) 



where the prime on the summation indic;atc that the term 
m = n is excluded from the sum. Next note that 



J q 

Substituting (A3) in (A2) gives 



477 



-^n I 



LHS 




{A3) 



If one assume that there is no correlation between the 
magnitude of (determined by the size of the contact) 
and its position, we can replace the individual current /„ 
in the double summation in (A2) by their mean and get 



LHS 



4n2 



-± — 

2^3 Iy... _ 



where n = N/Aq is the concentration of contact spots 
and N the total number of contact spots. 

In the same way as above one can simplify the expres- 
sion involving the normal stress (right hand side (RHS) 
of (Al)). We write 

n 

where /„ is the normal force acting in the contact n. 
Using this equation the RHS of (Al) becomes 



(AG) 



If one assume that there is no correlation between the 
magnitude of /„ and its position, we can replace the in- 
dividual current /„ in the double summation in (A6) by 
their mean and get 



RHS 



4n2 



-± — 

2_3 |y„_ _ 



^ I -^m -'^n I 



{A7) 



Thus, LHS « RHS and we have proved the (approxi- 
mate) equality (Al). 

Substituting (A5) in (20) gives 



1 1 



y— 



a TTKn N ^ Ixto — x„ 



(A8) 



which agree with the derivation of Greenwood 2l| . We 
refer to the article of Greenwood for an interesting dis- 
cussion about the contact resistance based on the (ap- 
proximate) expression (A8) for the contact resistance. 



This exact relation between the heat transfer coeffi- 
cient and the normal stiffness per unit area has already 
been derived by Barber [S^l using a someone different ap- 
proach. 



Appendix B 

The normal (interfacial) stress ^^(x) and the differ- 
ence in the surface displacement uoz(x) — uiz(x) at the 
interface can be considered to depend on the average in- 
terfacial separation u. The derivatives of these quantities 
with respect to u are denoted by cr^ and (j). In Appendix 
C we show that 



(q)-'5(q)-— AaUq). 



m 



Note that (15) and (Bl) are very similar. Thus, if we 
multiply both sides of (Bl) with AI and define M(j) = 
then (Bl) takes the form 



where 



^(q) = MS{q) - ^Aa'M 



2Mk 



E* 



Eq. (B2) is identical to (15) if we write 



or, equivalently, 



J^q) = /icr^,(q), 



J,(x)/Jo=ai(x)M 



iB2) 



(i?3) 



where Pq is the normal stiffness. We note that (B4) im- 
plies that the current density Jz (x) will be non- vanishing 
exactly where the normal stress (Tz(x) is non- vanishing, 
which must be obeyed in the present case, where all the 
heat current fiow through the area of real contact. How- 
ever, in order for Jz(x) to be proportional to crz(x) it is 
not enough that these functions obey similar (in the sense 
discussed above) differential equations, but both prob- 
lems must also involve similar boundary conditions. Now 
in the area of non-contact both and and hence cr^ 
must vanish. In the area of real contact the temperature 
field T is continuous so that ^ = T{x, -0)-T{x, +0) = 0, 
while the displacement field satisfies $ — uqz—uiz — h{x.) 
so that (since /i(x) is independent of u), (/) = in the area 
of real contact. Thus, both problems involves the same 
boundary conditions and Jz and a'^ must therefore be 
proportional to each other. 

Note that (B4) gives Jq = (J-Pq. Substituting (B3) in 
this equation and using the definition (16) gives an equa- 
tion of the form (3) with 

K dpo 



Appendix C 

In Ref. it was shown that the normal displacement 
uqz is related to the normal stress cr^ via 



wo^(q) = --^cFzio,), 

where = Eo/{l — Vq). In a similar way 

2 

uiziq) = ^crz(q). 



{CI) 



iC2) 



Let $ = uqz — uiz be the difference between the (inter- 
facial) surface displacement fields. Using (CI) and (C2) 
gives 



where 



$(q) = -^CTz(q) 
E*q 



11 1 



(C3) 



Note that the average of <&(x) is the average separation 
between the surfaces which we denote by u. Thus if 



we get 



az (x) = _po + Adz (x) 



$(q) -z25(q)- — Aaz(q) 



As the squeezing pressure po increases, the average sepa- 
ration u will decrease and we can consider po as a function 
of u. The quantity pg (u) is referred to as the normal stiff- 
ness per unit nominal contact area. Taking the derivative 
of (C4) with respect to u gives 



(q)='5(q)-^Aa;(q) 



(C5) 



where a'z is the derivative of az with respect to u and 
where = is the derivative of $ with respect to u. 

Appendix D 

Heat conduction result from the collisions between 
atoms as in fluids, or by free electron diffusion as pre- 
dominant in metals, or phonon diffusion as predominant 
in insulators. In liquids and gases, the molecules are usu- 
ally further apart than in solids, giving a lower chance of 
molecules colliding and passing on thermal energy. Met- 
als are usually the best conductors of thermal energy. 
This is due to the free-moving electrons which are able 



to transfer thermal energy rapidly through the metal. 
However, the difference in the thermal conductivity of 
metals and non-metals are usually not more than a fac- 
tor ^ 100. Typical values for the heat conductivity are 
K « 100 W/mK for metals, « 1 W/mK for insulators 
(e.g., metal oxides or polymers), ~ 0.1 W/mK for fluids 
(but for water n « 0.6 W/mK) and « 0.02 W/mK for 
gases at normal atmospheric pressure and room temper- 
ature. 

In contrast to thermal heat transfer, electric conduc- 
tion always involves the motion of charged particles (elec- 
trons or ions). For this reason the electric contact resis- 
tance is much more sensitive to oxide or contamination 
layers at the contacting interface then for the heat trans- 
fer. For the electric conduction the variation of the con- 
ductivity between good conductors (most metals), with 
the typical electric conductivity k' « 10^ (ilm)^^, and 
bad conductors such as silicon dioxide glass or (natu- 
ral) rubber where k' « 10^^^ {flm)~^, is huge. This 
makes the electrical contact resistance of metals sensi- 
tive to (nanometer) thin oxide or contamination layers. 
However, as pointed out in the Introduction, if there is a 
large number of small breaks in the film, the resistance 
may be almost as low as with no film. 

Appendix E 

Here we briefly summarize some results related to 
forced convective heat transferji^. When a fluid (e.g., 
air) flow around a solid object the tangential (and the 
normal) component of the fluid velocity usually vanish 
on the surface of the solid. This result in the formation 
of a thin boundary layer (thickness 5) at the surface of 
the solid where the fluid velocity rapidly increases from 
zero to some value which is of order the main stream 
velocity outside of the solid. If the temperature Ti at 
the solid surface is different from the fluid temperature 
Ifluid) the fluid temperature in the boundary layer will 
also change rapidly from Ti to Tfluid- Depending on the 
fluid flow velocity, the fluid viscosity and the dimension 
of the solid object the flow will be laminar or turbulent, 
and the heat transfer process is fundamentally different 
in these two limiting cases. In a typical case (for air) 
the thickness 5 ~ 1 mm and the heat transfer coefRcient 
a«K/(5« 10 W/m^K. 

Let us consider heat transfer from a rotating disk as a 
model for the heat transfer from a tire 43| . In this case 
it has been shown [3| that fully turbulent flow occur if 
the Reynolds number Re > 2.5 x 10^ where 



this case the heat transfer coefficient is given approxi- 
mately by 4J| : 



Re = 



vrR 



where R is the radius of the disk (or rather the distance 
from the center of the disk to some surface patch on the 
disk) , Lo the angular velocity and v the kinematic viscos- 
ity of air. In typical tire applications Re > 2.5 x 10^ so 
turbulent flow will prevail in most tire applications. In 



0.019- 



Sair / VrR 

R \ ly 



0.8 



As an example, at T = 300 K for air ly = 15.7x 10"*^ m^/s 
and Kair = 0.025 W/mK and assuming R — 0.3 m and 
vr — 30 m/s we get aair ~ 63 W/m^K. 
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